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Abstract 

Analyses of the inelastic a+^^C scattering at medium energies have indicated that 
the strength of the Hoyle state (the isoscalar 0^ excitation at 7.65 MeV in ^^C) 
seems to exhaust only 7 to 9% of the monopole energy weighted sum rule (EWSR), 
compared to about 15% of the EWSR extracted from inelastic electron scattering 
data. The full monopole transition strength predicted by realistic microscopic a- 
cluster models of the Hoyle state can be shown to exhaust up to 22% of the EWSR. 
To explore the missing monopole strength in the inelastic q+^'^C scattering, we 
have performed a fully microscopic folding model analysis of the inelastic a+^^C 
scattering at E'lab = 104 to 240 MeV using the 3-a resonating group wave function 
of the Hoyle state obtained by Kamimura, and a complex density-dependent M3Y 
interaction newly parametrized based on the Brueckner Hartree Fock results for 
nuclear matter. Our folding model analysis has shown consistently that the missing 
monopole strength of the Hoyle state is not associated with the uncertainties in the 
analysis of the a+^^C scattering, but is most likely due to the short lifetime and 
weakly bound structure of this state which significantly enhances absorption in the 
exit a+^^C*{0^) channel. 
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Given a vital role in the stellar synthesis of Carbon, the isoscalar O2 state at 
7.65 MeV in ^^C (known as the Hoyle state) has been studied over the years 
in numerous experiments. Although this state was clearly identified long ago 
in the inelastic a+^^C scattering at medium energies [T][2][S1H] and inelastic 
electron scattering [5j as an isoscalar EO excitation, our knowledge about its 
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unique structure is still far from complete [6] . Since the Hoyle state lies slightly 
above the a-decay threshold of 7.27 MeV, its wave function should have a dom- 
inant a cluster component. In fact, the Hoyle state has been well described 
in terms of three a clusters by the Resonating Group Method (RGM) some 
thirty years ago [TfH] . Very interesting is the condensate scenario suggested 
recently [9][T0] where the three a clusters were shown to condense into the 
lowest (s-state) of their potential, and thus, forming a Bose-Einstein conden- 
sate (BEG). A more complicated structure of the Hoyle state was shown by the 
Fermionic Molecular Dynamics (FMD) calculation [llj where the BEG wave 
function is mixed also with the molecular ®Be+a configuration. In general, to 
validate conclusion made in the structure calculation, the wave functions must 
be carefully tested in the study of nuclear reactions exciting the Hoyle state. 
In this aspect, the model wave functions of the Hoyle state Plllfl2j have been 
shown to give a reasonable description of the inelastic electron scattering data. 
The electric monopole transition moment 

oo 

M(EO, 0+ ^ 0+) = ev^ I p|;p4r)r^dr (1) 





predicted by the RGM |8], BEG [10] and EDM [11] calculations is around 6.62, 
6.45 and 6.53 e fm^, respectively, which are about 20% largeiQ than the ex- 
perimental moment M(i?0) ~ 5.37 ± 0.22 e fm^ deduced from the (e, e') data 
[5|. The inelastic electron scattering probes, however, only the charge transi- 
tion density and it is necessary to study also inelastic hadron scattering which 
probes the nuclear transition density. Such experiments have been done, e.g., 
for inelastic ^'*^He,^Li-|-^^G scattering p!3|14fl5j and the data analyses, either 
in the distorted wave Born approximation (DWBA) or coupled-channel (GG) 
formalism, were performed using the collective model form factor (EE) for the 
Hoyle state. These analyses indicated that the observed monopole transition 
strength of the Hoyle state exhausts about 7 to 9% of the isoscalar monopole 
energy weighted sum rule (EWSR). In particular, the DWBA analysis of the 
inelastic a+^^G data, measured recently at i^iab = 240 MeV with high preci- 
sion pJJ, has found that the Hoyle state exhausts 7.6 ± 0.9% of the EWSR. 
We recall that the sum rule fraction 5*0 of a monopole excitation is determined 
[IS] as 

S, = E^\M{ISO,Ot - 0+)|V C^^A^ll^] , (2) 



where E^, m and A are the excitation energy, nucleon mass and target mass 
number, respectively, < > is the mean square radius and M(JS'0, Oi" O2") 

^ a direct comparison of ([T]) with the observed M(i?0) can be made if the small 
difference between the charge- and proton transition densities is neglected. 
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is the isoscalar monopole transition moment which is determined by the same 
Eq. ([1]) but using the nuclear transition density instead of the charge transition 
density. If we assume M{ISO,Ot Ot) ~ 2M(E0,0^ 0^)/e and take 
< 5.78 fm^ (estimated using the ground state density of ^^C given by 

the RGM H]), then we obtain So ^ 22.8, 21.7 and 22.2% of the EWSR from the 
RGM [8], BEG [lO] and EDM [H] results, respectively, and Sq ^ 15.0±1.3% of 
the EWSR from the experimental monopole moment, in a perfect agreement 
with Sq value given in Ref. [5]. Given a good description of the (e, e') data 
by the cluster models pi5|lUllllj . we conclude that the monopole strength of 
the Hoyle state should exhaust 15 to 20% of the EWSR and a puzzle remains 
why at least half of the monopole strength of the Hoyle state is missing in the 
inelastic a+^^C scattering. 

To investigate the missing monopole strength of the Hoyle state, we have per- 
formed a consistent folding model analysis of the inelastic a+^^C scattering 
at -Eiab = 104 to 240 MeV using microscopic nuclear densities obtained from 
the 3-a RGM wave functions by Kamimura jB], and a new complex density- 
dependent M3Y interaction. Before discussing our results, we recall that the 
same RGM nuclear densities |8] were used in a number of folding model stud- 
ies (see, e.g., Refs. p!7|18|19ll20|21] ) of the inelastic a+^^C scattering. Except- 
ing the early work by Bauhoff [T7] where a density independent Gaussian 
has been used as aN interaction in the single-folding calculation, other stud- 
ies mentioned here have used the well-known DDM3Y \^2\ interaction in the 
double-folding calculation of the a+^^C potentials. Since the DDM3Y inter- 
action is real, the imaginary parts of both the optical potential (OP) and 
inelastic EE were chosen phenomenologically in these studies. Eor example, in 
the GC analysis of the a+^^G scattering by Ohkubo and Hirabayashi |18II19] 
the imaginary inelastic EE was neglected and parameters of the Woods-Saxon 
imaginary OP were adjusted separately for each exit channel to obtain a good 
GG fit to the measured cross sections. Although one could achieve a reason- 
able description of the inelastic a+^^G scattering data in such a GG analysis, 
an arbitrary choice of the imaginary potentials makes it difficult to estimate 
accurately the absolute EO transition strength. Up to now, the EO strength 
of the Hoyle state has been deduced from the DWBA or GG analyses of the 
inelastic ^'^He,^Li+^^G scattering [T3|14lll5j using the breathing mode (BM) 
model [16] for the nuclear transition density 

Po+->o+W = -/3[3po(r)+r^], (3) 

where po{r) is the ground-state (g.s.) density of ^^G. In this case, 5*0 = 
[P/Praax]"^, where Prnax IS the dcformatiou parameter required for a monopole 
excitation to exhaust 100% of the EWSR (see Eq. (3.5) in Ref. pj). The 
inelastic scattering EE is then obtained by (double) folding the transition 
density with the effective nucleon-nucleon (NN) interaction and projec- 
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tile g.s. density plpS] . In the a-nucleus scattering, one also uses a simpler 
approach to (single) fold the density ([3]) with an appropriate a-nucleon (aN) 
interaction [2B]. In general, the (complex) strength of the aN or A^A^ inter- 
action is first adjusted to the best optical model or CC description of elastic 
scattering and then is used without any further renormalization to calculate 
the inelastic FF. As a result, the only parameter in the inelastic channel is the 
deformation parameter (3 which must be determined from the best DWBA or 
CC fit to the inelastic scattering data. Since the dilute a-clustered structure 
of the Hoyle state is different from a monopole compression mode, the use 
of the BM density in the folding model calculation might be question- 
able. In addition, the uncertainties of the effective A^A^ or aN interaction and 
the optical potential (OP) in the entrance and exit channels could also cause 
some deviation in the description of inelastic cross section which results on 
the missing EO strength. 

We note that the RGM wave functions by Kamimura [S] were proven to give 
consistently a realistic description of the shell-like structure of the ground state 
Ot, first excited 2t (4.44 MeV) and SJ" (9.64 MeV) states, and the a-clustered 
structure of the 0^ (7.65 MeV) Hoyle state. For the latter, the RGM wave 
function has been shown [TUfT^ to be very close to the EEC wave function 
and both of them give, in fact, nearly the same description of the inelastic 
a+^^C scattering [20]. Thus, the RGM nuclear densities is a very good choice 
for the present folding model study. Concerning the effective A^A^ interac- 
tion, a density dependent version of the M3Y-Paris interaction (dubbed as 
CDM3Y6 interaction [21]) is the most often used in our folding calculations. 
The CDM3Y6 density dependent parameters were carefully adjusted in the 
Hartree-Fock (HF) calculation to reproduce the saturation properties of nu- 
clear matter [24j. However, like the DDM3Y version, the CDM3Y6 interaction 
is real and can be used to predict the real OP and inelastic FF only. To avoid 
a phenomenological choice of the imaginary potentials like those used in ear- 
lier folding model studies [T5|19II2U|[^ of the inelastic a+^^C scattering, we 
have constructed in the present work a new complex, density dependent M3Y- 
Paris interaction. The parameters of the complex density dependence were 
calibrated against the Brueckner Hartree-Fock (BHF) results [27] for the nu- 
cleon OP in nuclear matter by Jeukenne, Lejeune and Mahaux (JLM) at each 
considered energy. Namely, the isoscalar complex nucleon OP in nuclear mat- 
ter is determined from the HF matrix elements of the new density dependent 
A^A^ interaction as 

Uo{E,p) = J2 [< kj\uu{E,p)\kj > + < kj\uEx{E,p)\jk >]. (4) 

j<kF 

Here kp = [I.Stt^p]^/^ and k is the momentum of the incident nucleon which 
must be determined self-consistently from Uo ask = J2m[E — Re Uq{E, p)\/K^. 
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We have used in Eq. (jl]) two different CDM3Y functionals [23] to construct 
separately the real and imaginary parts of the density dependence, so that the 
real and imaginary parts of the effective A^A^ interaction 'Ud(ex) are determined 
as 



Re[MD(EX)] = Fy{E,p)vY)(^EX){s), Im[MD(EX)] = F^^{E , p)vr,{EX)is) . (5) 



The radial direct and exchange interactions fD(EX)(s) were kept unchanged, 
as derived from the M3Y-Paris interaction [2S|, in terms of three Yukawas 
(see, e.g., Ref. [25||29] ). The parameters of the complex density dependence 
Fv(w)(-E', p) were adjusted iteratively until Uq{E,p) agrees closely with the 
tabulated JLM results at each energy [27] as shown in Fig. [H More details on 




Fig. 1. Isoscalar optical potential Uq{E,p) of nucleoli incident on nuclear matter at 
43 and 60 MeV given by HF calculation ^ using the new CDJLM interaction (solid 
curves). The points are the microscopic BHF results made available up to nuclear 
matter density p = 0.2 fm~'^ by the JLM group ^27j. This CDJLM interaction has 
been used in the present folding model study of the a+^^C scattering at ^lab = 172.5 
and 240 MeV. 

the new density dependent interaction, dubbed hereafter as CDJLM interac- 
tion, and the explicit parameters of Fv(w)(-E', p) at different energies will be 
presented elsewhere. The CDJLM interaction was used to calculate the com- 
plex OP and inelastic FF for the microscopic DWBA and CC analyses of the 
elastic and inelastic a+^^C scattering data at i?iab = 104 [Tp] . 139 [3], 172.5 
[4j, and 240 MeV [15j. The generahzed double- folding method [25] was used 
to calculate the complex a-nucleus potential as the following HF-type matrix 
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element of the CDJLM interaction ([5]) 



[< U>d|u > + < ifluExlji >], 



(6) 



ieay'eAj'eA* 



where A and A* are states of the target in the entrance- and exit channel of the 
a-nucleus scattering, respectively. Thus, Eq. IQ gives the elastic (diagonal) 
OP if A* = A and inelastic (transition) FF if otherwise. A more accurate local 
density approximation suggested in Ref. [30] has been used for the exchange 
term in Eq. (Q. The dynamic change in the density dependence Fv(w)(p) 
caused by the inelastic a+^^C scattering [25|3T] is taken into account prop- 
erly. All the DWBA and CC calculations have been performed using the code 
EC1S97 written by Raynal 



Fig. 2. Elastic a+^^C scattering data measured at Ei^^y, = 104 pj2j, 139 ^, 172.5 [4J 
and 240 MeV [15j in comparison with the OM and CC results given by the complex 
folded OP. The elastic a scattering on the Hoyle state is predicted by the OM 
calculation using the diagonal a+^^C*{0^) folded OP with Ni chosen (see Table [I]) 
to reproduce the measured inelastic a+^'^C*(0^) scattering cross section in the CC 
calculation. 

To fine tune the strength of complex CDJLM interaction ([5]) which has been 
obtained in the nuclear matter limit (jl]), the optical model (OM) analysis of 
the measured elastic a+^^C scattering data was done at each energy using the 
complex OP determined from the elastic folded potential ([6]) as 




Elastic scattering 
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Uo{R) = NnRe[U,+^,+ iR)] + iNjlm[U,+^,+ {R)] + VdR), 



(7) 



6 



where Vc{R) is the elastic Coulomb potential taken, for simplicity, as that 
between a point charge and a uniform charge distribution of radius Rq = 
1.3 X (4^/'^ + 12^/^) fm. The renormalization coefficients iV^ and A^i of the real 
and imaginary elastic folded potentials were adjusted first to the best OM 
fit of the elastic scattering data and then used further to scale the real and 
imaginary inelastic folded FF for the DWBA calculation, a standard method 
used sofar in the folding + DWBA analyses of inelastic a-nucleus scattering 
|15|17|25ll26] . Since N-^ and Ni are an approximate way to take into account 
the higher-order (dynamic polarization) contributions to the microscopic OP 
they need to be readjusted again (to fit the elastic data) in the CC 
calculation, to account for the remaining nonelastic channels which were not 
included into the CC scheme. For consistency, these new Nn and Ni factors 
are also used to scale the complex inelastic folded FF for the CC calculation. 
The OM and CC descriptions of the elastic a+^^C scattering at the considered 
energies are shown in Fig. [2] and the corresponding OP parameters are given 
in Tabled! One can see from Tabled] that the A^r coefficient given by the best 
OM fit of the elastic scattering data is slightly above unity which indicates 
that the new CDJLM interaction is quite a realistic choice for the double- 
folding calculation of the a-nucleus potential. The best fit A^i of about 1.3 to 
1.5 given by the OM calculation is not unexpected because the imaginary part 
of the CDJLM interaction is based on the BHF results for nuclear matter and 
gives, therefore, only a "volume" absorption. As a result, the imaginary folded 
OP cannot properly account for the surface absorption caused by inelastic 
scattering to the low-lying collective excitations and transfer reactions, and 
the OM fit to the elastic data naturally requires an enhanced A^i coefficient. 
The best-fit folded optical potentials give volume integrals Jr{i) and total 
reaction cross sections ctr very close to those obtained earlier in the folding 
model analysis using a phenomenological imaginary OP [30] or in the model- 
independent Fourier-Bessel analyses of the elastic a+^^C scattering data [33j- 
Thus, the elastic distorted waves given by the CDJLM folded OP should be 
quite accurate for the DWBA analysis of the inelastic a+^^C scattering. 

The DWBA results obtained with the complex folded OP and inelastic folded 
FF are compared with the measured cross section for inelastic a+^^C scatter- 
ing to the Hoyle state in Figs. [3] and IH We found that the calculated DWBA 
cross sections systematically overestimate data at all energies if the inelas- 
tic folded FF is obtained with the {unrenormalized) RGM transition density 
[8]. Given the accurate choice of the density dependent A^A^ interaction and 
nuclear densities, this discrepancy is clearly not associated with usual uncer- 
tainties of the folding model analysis of a-nucleus scattering [2n]- Since the 
RGM transition density has been proven to reproduce nicely the (e, e') data 
[81IT2] , the folding + DWBA results shown in Figs. [3] and H] might indicate a 
"suppression" of the monopole EO strength occurred in the inelastic a+^'^C 
scattering. To further probe this effect, we have used a realistic Fermi distri- 
bution for the g.s. density of ^^C [20] to generate the BM transition density 
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Table 1 

Renormalization coefficients A'j^(j) of the complex folded OP and inelastic FF used 
in the DWBA and CC analyses of a+^^C scattering at ^lab = 104, 139, 172.5, 
and 240 MeV. </r(i) are volume integrals of the real and imaginary OP, and cjr is 
the corresponding total reaction cross section. CCen are the OP parameters in the 
entrance a+^^Cg.s. channel of the CC calculation. CCex are those of the OP in the 
exit a+^^C*(Oj) channel, where A^i was adjusted to reproduce the inelastic 0^ data 
using the original ROM transition density [8j. 



^lab (MeV) 


Analysis 




Jr (MeV fm3) 




Ji (MeV fm^) 


cjR (mb) 


104 


DWBA 


0.975 


312.3 


1.295 


104.3 


799.5 




CCcn 


1.015 


325.1 


1.066 


81.76 


794.1 




CCex 


1.015 


325.1 


2.500 


201.4 


1349 


139 


DWBA 


1.025 


287.4 


1.374 


120.2 


770.8 




CCcn 


1.049 


294.2 


1.064 


93.06 


739.7 




CCex 


1.049 


294.2 


2.500 


218.6 


1287 


172.5 


DWBA 


1.156 


285.4 


1.506 


117.5 


736.7 




CCcn 


1.165 


287.6 


1.264 


98.64 


716.9 




CCex 


1.165 


287.6 


3.000 


234.1 


1288 


240 


DWBA 


1.127 


252.8 


1.340 


111.0 


659.1 




CCcn 


1.145 


256.9 


1.241 


102.8 


656.4 




CCex 


1.145 


256.9 


3.400 


281.6 


899.1 



for the DWBA analysis of inelastic scattering to the Hoyle state. We have fixed 
the monopole deformation parameter j3 so that gives exactly the same EO 
transition strength as that given by the RGM transition density and, hence, 
exhausts 22.8% of the isoscalar monopole EWSR as discussed above. As can 
be seen from Figs. [3] and HJ the DWBA cross sections given by the inelastic 
FF folded with the BM transition density for the Hoyle state are very close to 
those given by the RGM transition density. To match the calculated DWBA 
cross section to the data points, we need to scale both transition densities 
by a factor of 0.55 which reduces the sum rule strength of the Hoyle state 
to 5*0 ~ 6.9% of the EWSR. Such a small sum rule strength is significantly 
below the empirical range of 15 - 20% of the EWSR discussed above. Given 
an accurate energy dependence of the CDJLM interaction, with its density 
dependence carefully calibrated against the BHF results at each energy, the 
RGM and BM transition densities scaled by the same factor of 0.55 give rea- 
sonable DWBA descriptions of the data at four considered energies. The sum 
rule strength of the Hoyle state found in our DWBA analysis is comparable 
to that (7.6 ± 0.9%) found in the folding + DWBA analysis of the inelastic 
a+^^C scattering at 240 MeV [U], where a Gaussian has been used for the 
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Fig. 3. Inelastic a+^^C scattering data at ^lab = 104 [1|2] and 139 [3] MeV for 
the 0^ state of in comparison with the DWBA and CC results given by the 
complex folded OP and inelastic FF obtained with two choices of the transition 
density for the O2 state. A good DWBA description of the data is reached only if 
the the transition densities are scaled by a factor of 0.55 which reduces EO sum rule 
strength to Sq ~ 6.9% of the monopole EWSR. If So is kept at 22.8% of the EWSR, 
as given by the RGM result [8], the data can be reproduced in the CC calculation 
only by using a more absorbing OP for the a+^'^C*(0^) channel. See more details 
in text. 

aN interaction in the folding calculation of the OP and inelastic FF. 

A widely used prescription for the DWBA calculation of inelastic hadron scat- 
tering is to assume the same complex OP for the entrance and exit channels 
[T6] and the only difference in the elastic waves is a slight shift in kinematics 
caused by a lower energy in the exit channel. For such a subtle nuclear state 
like the Hoyle state which has unusually dilute and extended structure [S] , the 
diagonal potential Uq+^q+ should be quite different from Uq+^q+ which was 
used as the OP in both the entrance and exit channels of the DWBA calcu- 
lation presented above. Actually, the importance of using explicitly different 
OP for different exit channels in the CC description of a+^^C scattering at 
104 and 139 MeV was pointed out earlier by Bauhoff [T7]. Thus, one must go 
beyond DWBA and do the CC analysis of inelastic a+^^C scattering to the 
Hoyle state before making any conclusion on the missing monopole strength. 
Although the E2 transition between 0^ and 2f states is quite weak and 0^ 
state does not belong to a two-phonon (2^ 2^) band, given the importance 
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Fig. 4. The same as Fig. [3] but for the data measured at E'lab = 172.5 [4J and 240 
MeV. 

of 0^ —>■ 2f transition in the stellar Carbon production, we have considered 
following coupling scheme in our CC calculation 

0+^2+^0^^ 0^ (8) 



In the CC scheme ([H]), the two-step excitation {Oi ^ 2^ ^ 0^) of the Hoyle 
state is treated in equal footing with the direct excitation {Of — > Thus, the 
CC description of the inelastic a+^^C scattering to state should be more 
realistic compared to the DWBA calculation which takes into account the 
direct excitation only [35]. For the inputs of the CC calculation, in addition 
to three inelastic FF's for transitions between Of, 2i and states, three 
complex OP's have been calculated using the diagonal nuclear densities of 
these states given by the ROM calculation [8] and CDJLM interaction ([5]). 
For consistency, the A^r(i) factors have been readjusted again to obtain a good 
CC description of the elastic scattering data, as shown in Fig. [2l While the 
best-CC-fit iVpi, factor (see Table [1]) remains more or less the same, A^i factor 
becomes significantly smaller and closer to unity, especially at the two low 
energies. This result shows explicitly how the inelastic channels in the CC 
scheme ([8]) contribute to the enhanced A^i factor found in the OM analysis 
of elastic scattering. With the increasing energy, more nonelastic channels 
are open and the best-CC-fit Ni factor remains around 1.2 to 1.3 at 172.5 
and 240 MeV. All the CC results are shown in the right panels of Figs. [3] 
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and m and it can be seen that the CC cross sections strongly overestimate the 
inelastic data for 0^ state, in the same way as found in the DWBA calculation, 
when the inelastic FF is obtained with the original RGM transition density 
[S]. The two-step excitation {Of ^ 2f ^ Of) was found to be negligible 
at all considered energies and the CC cross sections shown in Figs. [3] and 
m are practically given by the direct (one-step) excitation of state only. 
We note that the original RGM density for the 2^ —>■ Of transition could 
not reproduce the observed B{E2) transition probability and we have scaled 
this transition density, as recommended in Ref. [8], by a factor of 1.53 to 
reproduce the experimental transition rate B{E2, 2f Of) ~ 2.6 ±0.8 e^fm^ 
[36] . Despite such an enhancement, the folded 2f Of transition FF gives 
still a very weak contribution to the CC cross section for the excitation. 
An agreement with data can be reached only if the EO transition density is 
reduced by about half as found in the DWBA analysis. It is clear now that 
the missing monopole strength in the inelastic a+^^C scattering to the Hoyle 
state is not associated with uncertainties of the DWBA or CC analyses, and 
there should be some physics effect that damps the EO transition strength. 
We further mention a recent folding + CC calculation of inelastic a+^^C 
scattering at Ei^h = 172.5, 240 and 386 MeV (using the same RGM nuclear 
densities and DDM3Y interaction) by Takashina [37], where coupling to all 
the low-lying excited states of ^^C as well as 3-a breakup channel was taken 
into account, but the CC cross section for the 02" excitation remains lying 
significantly higher than the data points, in about the same manner as found 
in our CC analysis. 

Such a "damping" of the transition strength in inelastic nucleus-nucleus scat- 
tering has also been observed in our recent folding model study [38] of inelastic 
^^O+^^O scattering to the 2f and 3j~ states of ^^O. Given the electric transition 
strengths of these states of ^^O well determined from the (e, e') data (like the 
Hoyle state considered here), the DWBA or CC description of high-precision 
inelastic ^^O+^^O scattering data, which cover a wide angular range and 6 or- 
ders of the cross-section magnitude, was possible only if the absorption in the 
exit channel is significantly increased [3^. Such an enhanced absorption was 
also found for the exit channel of one-neutron transfer 0l7ol5 reaction to the 
3/2~ excited state of ^^O [39]. The enhanced absorption in the exit channel 
is a direct consequence of the suppression of nuclear refraction in nonelastic 
channels [29]. Namely, a short-lived (excited) and loosely bound cluster like 
^^□2+ or ^^03^2- has a shorter mean free path in the nuclear medium which 
implies a weaker refraction or stronger absorption [40|l4T] in the exit chan- 
nel compared to the entrance channel where both nuclei are in their (stable) 
ground states (for detailed discussion see Refs. [29|38f39j ). In a similar sce- 
nario, the 02" state of ^^C is very weakly bound and particle unstable (with 
a mean lifetime r ~ 10~^^ s) and it is natural to expect that the OP in 
the a+^^C*(02') exit channel is more absorbing than the OP in the entrance 
a+^^C channel. Since the elastic a scattering on the Hoyle state cannot be 
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measured, we have adjusted the renormahzation factor A^i of the imaginary 
OP of the q;+^^C*(0^) channel to achieve a reasonable CC description of the 
inelastic data for 0^ state, while keeping A^r factor and complex strength of 
the transition FF unchanged. From Figs. [3] and H] one can see that a very good 
CC description of the inelastic 0^ data can be reached by using A^i ^ 2.5 — 3.4 
for the imaginary OP in the a+^^C*(02') channel (see Table [1]). Such a strong 
absorptive OP predicts a very different elastic a+^^C*(02') scattering cross 
section compared to the measured elastic a+^^Cg.s. scattering (see Fig. [2]). 




®c...(deg) 

Fig. 5. Inelastic q+^^C scattering data measured at -Eiab = 104 |l|2j . 139 [3], 172.5 
[3] and 240 [15] MeV for the 2f excitation at 4.44 MeV in ^^C in comparison with 
the DWBA and CC results obtained with the complex folded OP and inelastic 
folded FF. See more details in text. 

To illustrate whether the enhanced absorption is also required for some other 
excited states of ^^C or is it a unique effect associated with fragile structure 
of the Hoyle state, we have done the same folding + DWBA (CC) analyses 
of inelastic a+^'^^C scattering to the 2^ state at 4.44 MeV using the ROM 
density and the results are plotted in Fig. [51 Like the Hoyle state, the RGM 
transition density of the 2f state has also been used to successfully reproduce 
the (e, e') data [8j. The electric transition strength predicted by the RGM 
transition density is B[E2, Of — > 2f) = 46.5 e^fm^ which agrees fairly with 
the measured B{E2) transition rate of 40 ± 4 e^fm^ [32112] • We note that 
2f state is a very strong excitation and the measured inelastic cross section 
becomes even comparable to the elastic scattering cross section at the medium 
angles where the nuclear scattering dominates (see Figs. [2] and [5]). That is the 
reason why the coupling strength by the 2^ excitation is strongest in our CC 
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scheme (|8]). From Fig. [5]one can see that a consistent folding model description 
of the inelastic a+^^C scattering to the 2i state at a-energy up to 172.5 MeV 
can be reached only in the CC calculation. If one stays within the DWBA 
framework, using the same inelastic folded FF, then the total strength of the 
E2 transition density needs to be enlarged to reproduce the 2f data, but this 
would lead to a larger discrepancy between the calculated and experimental 
B{E2,0f —>■ 2f) values. It can also be seen from Table [1] that the coupling 
strength by the 2i excitation exhausts nearly all the surface enhancement 
of the imaginary OP at 104 and 139 MeV, and the best-CC-fit Ni factor is 
around unity compared to that of about 1.3 - 1.4 as given by the OM analysis. 
Finally, the measured inelastic a+^^C scattering cross section for the 2^ state 
is nicely reproduced by the CC calculation at all energies without the need 
to enhance the imaginary OP in the a+^'^C*{2i) channel, i.e., CCen set of 
the iVR(i) factors in Table [1] was also used for the exit a+^^C*(2^) channel. 
This indicates that the enhanced absorption found above for the a+^^C*(0^) 
channel seems to be associated particularly with the fragile structure and 
short lifetime of the Hoyle state. We recall that the 2f state of ^^C is more 
robust and particle stable, and it decays via 7-emission only. With a mean 
lifetime r ~ 6 x 10~^^ s [42j the 2f state lives about 600 times longer than 
the Hoyle state and 10 times longer than the 2f state of ^^O for which the 
absorption enhancement has been found in the folding model study of inelastic 
160+160 scattering p]. Although our folding model analysis was done using 
a new complex density dependent interaction, the absorption enhancement 
found for the exit a+^^C*(0^) channel is not a feature associated with the 
new interaction. In fact, such a large absorption has already been established 
in the semi-microscopic CC analysis of the same a+^^C data by Ohkubo and 
Hirabayashi [19]. These authors have used phenomenological Woods-Saxon 
(WS) potentials of the same geometry for the imaginary OP in all exit channels 
of the inelastic a+^^C scattering and adjusted the WS depth each case to 
reproduce the inelastic data. From Table I of Ref. [19] one can find that the 
WS imaginary OP for the a+^^C*(0^) channel is nearly 4 times deeper than 
that for the entrance a+^^Cg.s. channel, while the best-fit WS depth found 
for the a-t-^^C*(2^) channel is only slightly larger than that of the a-|-^^Cg,s. 
channel. For the 3i state of ^^C the CC calculation of Ref. \i9l has shown also 
a sizable enhancement of absorption, in an agreement with the results of our 
CC analysis using the CDJLM interaction which will be presented elsewhere. 
We note further that small fraction (5*0 ~ 9%) of the monopole EWSR found 
in the DWBA analyses of inelastic ^He,^Li+^^C scattering to the Hoyle state 
[T3|IT4] could well be associated with the same absorption effect as that found 
in the present work. 

Our results of the CC description of inelastic a+^^C scattering to the Hoyle 
state as well as those of inelastic scattering and one-neutron transfer reac- 
tion measured with the ^^0-1-^^0 system [381139] show clearly that there is 
a correlation between the weak binding and/or short lifetime of the excited 
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target-like cluster and the absorption in the exit channel of a quasi-elastic 
scattering reaction. Although it is still difficult to establish a systematics for 
this "absorption" phenomenon, the results obtained sofar indicate the need 
to have a realistic choice for the OP not only in the entrance but also in the 
exit channel of the inelastic nucleus-nucleus scattering. The standard method 
of choosing the same complex OP in both the entrance and exit channels of 
nucleus-nucleus scattering might lead to a large uncertainty in the deduced 
transition strength if one simply scales the inelastic FF to match the DWBA 
or CC results to the measured angular distributions. In particular, this effect 
must be accurately taken into account in the study of the nucleus-nucleus 
scattering measured with unstable beams where nuclei in the entrance and 
exit channels are very differently bound. 



EO transition density for Hoyle state 




Fig. 6. OJ" — > O2 transition densities for the Hoyle state given by the RGM calcu- 
lation of Kamimura [8j and BM model The deformation parameter P of the 
BM density was adjusted to reproduce the same monopole transition strength as 
that given by the RGM transition density which exhausts 22.8% of the isoscalar 
monopole EWSR. 

We discuss now briefly the choice for the EO transition density of the Hoyle 
state. It is a common belief that the breathing mode (BM) model ([3]) is 
more appropriate for the transition density of a giant monopole resonance 
(in medium and heavy nuclei) that exhausts a large fraction of the monopole 
EWSR [23] ■ Although the BM density was scaled to reproduce the monopole 
transition moment given by the RGM density [S] before being used in our 
folding calculation, its radial dependence (see Fig. [6]) should not contain any 
feature of the a cluster or BEG structure of the Hoyle state. Therefore, it is 
quite surprising to see that the BM transition density gives nearly the same 
DWBA and GG descriptions of the inelastic data as those given by the 3-a 
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RGM transition density (see Figs. [3] and H]). One can see in Fig. |6] that the 
two densities have the same asymptotic tail at r > 6 fm while the node and 
two extrema of the BM density are shifted inward to smaller radii compared 
to the RGM density. However, as shown in Figs. [3] and HJ such a difference in 
the two densities gives only a minor difference in the calculated inelastic scat- 
tering cross sections. Consequently, the success of the BEC or RGM nuclear 
densities in the CC description of inelastic a+^^C scattering to the 0^ state 
reported sofar [T8|19II20||2T] is an important evidence, but not the unambigu- 
ous proof for the a-condensate structure of the Hoyle state. In fact, a closer 
inspection of the structure models by Chernykh et al. [TT] has revealed that 
the structure of the Hoyle state could be somewhat more complicated than 
just a condensate of three a particles. 

In conclusion, a realistic (complex) density dependence was introduced into 
the M3Y-Paris interaction, based on the Brueckner Hartree Fock calculation 
of nuclear matter, for the folding model study of the a+^^C scattering at 
-^'lab = 104 to 240 MeV. Given an accurate estimation of the bare a+^'^C 
optical potential, our folding model analysis has shown consistently that there 
should be an enhancement of absorption in the exit a-|-^^C*(02') channel due 
to the short lifetime and weakly bound structure of the Hoyle state, which 
accounts for the "missing" monopole strength of the Hoyle state observed 
earlier in the DWBA analysis of inelastic a+^^C scattering. 



Acknowledgments 

The authors are grateful to Prof. W. von Oertzen, Prof. S. Ohkubo, Dr. M. 
Takashina and Dr. X. Chen for their helpful communications. We also thank 
Prof. M. Kamimura for providing us with the revised parametrization of the 
RGM densities from Ref. 0. 



References 

[1] G. Hauser, R. Lohken, H. Rebel, G. Schatz, G.W. Schweimer, and J. Specht, 
Nucl. Phys. A 128 (1969) 81. 

[2] J. Specht, G.W. Schweimer, H. Rebel, G. Schatz, R. Lohken, G. Hauser, Nucl. 
Phys. A 171 (1971) 65. 

[3] S.M. Smith, G. Tibell, A.A. Cowley, D.A. Goldberg, H.G. Pugh, W. Reichart, 
N.S. Wall, Nucl. Phys. A 207 (1973) 273. 



15 



[4] S. Wiktor, C. Mayer-Boricke, A. Kiss, M. Roggc, P. Turek, and H. Dabrowski, 
Acta Phys. Pol. B 12 (1981) 491; A. Kiss, C. Mayer-Boricke, M. Rogge, P. 
Turek, and S. Wiktor, J. Phys. G 13 (1987) 1067. 

[5] P. Strehl, Z. Phys. 234 (1970) 416. 

[6] M. Freer, J. Phys. G 34 (2007) 789; M. Freer, to appear in Proceedings of Intern. 
Symp. on Physics of Unstable Nuclei (ISPUN07) (World Scientific: Singapore). 

[7] E. Uegaki, S. Okabe, Y. Abe, H. Tanaka, Prog. Theor. Phys. 57 (1977) 1262. 

[8] M. Kamimura, Nucl. Phys. A 351 (1981) 456; M. Kamimura, private 
communication . 

[9] A. Tohsaki, H. Horiuchi, P. Schuck, G. Ropke, Phys. Rev. Lett. 87 (2001) 
192501. 

[10] Y. Funaki, A. Tohsaki, H. Horiuchi, P. Schuck, G. Ropke, Phys. Rev. C 67 
(2003) 051306(R). 

[11] M. Chernykh, H. Feldmeier, T. Neff, P. von Neumann-Gosel, A. Richter, Phys. 
Rev. Lett. 98 (2007) 032501. 

[12] Y. Funaki, A. Tohsaki, H. Horiuchi, P. Schuck, G. Ropke, Eur. Phys. J. A 28 
(2006) 259. 

[13] D. Lebrun, M. Buenerd, P. Martin, P. de Saintignon, G. Perrin, Phys. Lett. 97 
B (1980) 358. 

[14] W. Eyrich, A. Hofmann, A. Lehmann, B. Miihldorfer, H. Schlosser, H. Wirth, 
H. J. Gils, H. Rebel, S. Zagromski, Phys. Rev. C 36 (1987) 416. 

[15] B. John, Y. Tokimoto, Y.W. Lui, H.L. Clark, X. Chen, D.H. Youngblood, Phys. 
Rev. C 68 (2003) 014305. 

[16] G.R. Satchler, Direct Nuclear Reactions (Clarendon Press: Oxford, 1983). 

[17] W. Bauhoff, Phys. Lett. 139 B (1984) 223. 

[18] Y. Hirabayashi, S. Ohkubo, Phys. At. Nuclei 65 (2002) 683. 

[19] S. Ohkubo, Y. Hirabayashi, Phys. Rev. C 70 (2004) 041602(R). 

[20] M. Takashina, Y. Sakuragi, Phys. Rev. C 74 (2006) 054606. 

[21] S. Ohkubo, Y. Hirabayashi, Phys. Rev. C 75 (2007) 044609. 

[22] A.M. Kobos, B.A. Brown, P.E. Hodgson, G.R. Satchler, A. Budzanowski, Nucl. 
Phys. A 384 (1982) 65. 

[23] D.J. Horen, J.R. Beene, G.R. Satchler, Phys. Rev. C 52 (1995) 1554. 

[24] D.T. Khoa, G.R. Satchler, W. von Oertzen, Phys. Rev. C 56 (1997) 954. 

[25] D.T. Khoa, G.R. Satchler, Nucl. Phys. A 668 (2000) 3. 



16 



[26] G.R. Satchler, D.T. Khoa, Phys. Rev. C 55 (1997) 285. 

[27] J.R Jeukenne, A. Lejeune, C. Mahaux, Phys. Rev. C 16 (1977) 80. 

[28] N. Anantaraman, H. Toki, and G.F. Bertsch, Nucl. Phys. A 398 (1983) 269. 

[29] D.T. Khoa, W. von Oertzen, H.G. Bohlen, S. Ohkubo, J. Phys. G 34 (2007) 
Rill. 

[30] D.T. Khoa, Phys. Rev. C 63 (2001) 034007. 

[31] D.K. Srivastava, H. Rebel, J. Phys. G 10 (1984) L127. 

[32] J. Raynal, Computing as a Language of Physics (IAEA, Vienna, 1972) p.75; J. 
Raynal, coupled-channel code ECIS97 (unpublished). 

[33] C.J. Batty, E. Friedman, H.J. Gils, H. Rebel, Adv. Nucl. Phys. 19 (1989) 1. 

[34] H. Feshbach, Theoretical Nuclear Physics Volume II (Wiley-Inter science: New 
York, 1992). 

[35] W. von Oertzen, private communication. 

[36] P.M. Endt, At. Data and Nucl. Data Tables 23 (1979) 3; F. Ajenberg-Selove, 
Nucl. Phys. A 506 (1990) 1. 

[37] M. Takashina, talk presented at DREB07 Workshop, June 2007. 

[38] D.T. Khoa, H.G. Bohlen, W. von Oertzen, G. Bartnitzky, A. Blazevic, F. 
Nuoffer, B. Gebauer, W. Mittig, P. Roussel-Chomaz, Nucl. Phys. A 759 (2005) 
3. 

[39] H.G. Bohlen, D.T. Khoa, W. von Oertzen, B. Gebauer, F. Nuoffer, G. 
Bartnitzky, A. Blazevic, W. Mittig, P. Roussel-Chomaz, Nucl. Phys. A 703 
(2002) 573. 

[40] J.P. Jeukenne, A. Lejeune, C. Mahaux, Phys. Rep. 25 (1976) 83; Phys. Rev. C 
16 (1977) 80. 

[41] J.W. Negele, K. Yazaki, Phys. Rev. Lett. 47 (1981) 71. 

[42] S. Raman, C.W. Nestor Jr., P. Tikkanen, At. Data and Nucl. Data Tables 78 
(2001) 1. 



17 



